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Abstract Mathematical analysis of the Anderson localization has been 
facilitated by the use of suitable fractional moments of the Green function. 
Related methods permit now a readily accessible derivation of a number 
of physical manifestations of localization, in regimes of strong disorder, 
extreme energies, or weak disorder away from the unperturbed spectrum. 
The present work establishes on this basis exponential decay for the mod- 
ulus of the two-point function, at all temperatures as well as in the ground 
state, for a Fermi gas within the one-particle approximation. Different im- 
plications, in particular for the Integral Quantum Hall Effect, are reviewed. 
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1. Introduction 



l.a The localization condition 

This is a report on recent progress in the mathematical analysis of Anderson local- 
ization. The simplifications which have been made in its derivation permit now to access 
a number of interesting properties of systems with disorder, by methods which are both 
mathematically rigorous and not excessively complicated. The report includes some new 
technical statements, but we recall also a number of previously known results, derived 
by various other authors, in order to present a more complete picture of the physically 
motivated questions which can be addressed by related mathematical methods. 

Anderson localization was first discussed in the context of the conduction proper- 
ties of metals ([1, 2]), but the mechanism is of relevance in a variety of other situations 
(e.g., [3]). The basic phenomenon is that disorder can cause localization of electron 
states (or normal modes — in other systems) and thereby affect properties like: time 
evolution (non-spreading of wave packets), conductivity (in response to electric field), 
Hall currents (in the presence of both magnetic and electric field) , and statistics of the 
spacing between nearby energy levels. 

In the electron gas approximation the system of electrons in a crystal is modeled 
by a gas of Fermions moving on a lattice. We focus here on systems with homogeneous 
disorder, which otherwise are periodic or translation invariant, at least up to gauge 
transformations. The excitations of the system are described by an effective one-body 
Hamiltonian, which consists of a short-range hopping term and a local potential. The 
one-particle Hamiltonian is a self-adjoint operator with matrix elements of the form 

H = K Xjy + UF r + \V x , (1.1) 

acting in the Hilbert space £ 2 (Z d ), where K x , y is a short-range hopping term, U^ er a 
periodic potential, and AV" a random potential expressing the disorder (impurities) with 
a tunable strength parameter A. 

We shall not discuss here the validity of the one-particle approximation, or that of 
the linear response theory. Instead, we focused on the analysis within such frameworks. 
In particular, we shall demonstrate how resolvent estimates can be used to address a 
number of physically motivated questions. 

For K x y we consider the following two cases: 

No magnetic field: K XtV depend only on the difference (x — y). 

Constant magnetic field: There is some ambiguity in the definition of the magnetic flux, 
since flux differences of h/e can be induced, or compensated for, by gauge transforma- 
tions. For concreteness sake, let us restrict to the operators of the form: 

K x>y = e~ iA ™ 6\ x -y\ tl (1.2) 
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with a phase A x , y which is an antisymmetric function of the oriented bonds, b = {x, y}. 
(A Xt y can be viewed as the line integral of the "vector potential" x(—e/h) along the 
direct path from x to y). The magnetic flux through a plaquette P is taken to be 

B P = -- £ Arg(^), (1.3) 
e bedP 

with the argument function interpreted through its principal branch, i.e., — n < Argz < 
7r. At non-zero field, translation invariance is possible only in the sense of magnetic 
translations, which combine shifts with gauge transformations, i.e. are unitaries of the 
form 

U(a)\x) = e-' lLpa{x) \x + a) . (1.4) 

In such cases, K = U (a)KU (a)* implies ordinary translation invariance for gauge invari- 
ant quantities, such as \K XjV \ and \ (x\(K — E — ir])~ 1 \y)\. (The fact that the composition 
law for the magnetic shifts provides only a projective representation of the translation 
group does not affect our analysis.) 

The potential V is realized as a collection of independent identically distributed 
random variables V x , whose probability distribution may be of the form r(v)dv with 
r(v) a bounded probability-density function. (These conditions may be relaxed: the 
results described below are valid also for a broad class of correlated randomness, more 
singular probability distributions for V, and Hamiltonians with off diagonal disorder, 
i.e., randomness in K xy .) 

Of central importance in the analysis is the Green function, i.e., the kernel of the 
resolvent operator: 

G(x,y;E + ir ] ) = (x\ H _ 1 E _^ \y) . (1.5) 

The behavior of this function at rj = 0+ reveals a great deal about the spectral prop- 
erties of the Hamiltonian (e.g., discrete versus continuous spectrum), the nature of its 
eigenfunctions (localized or extended), and the response of the system (e.g., to electric 
fields) at the linear-response level. 

A technically convenient signature of localization is a bound on the fractional 
moments of G(0, x; E). The explicit condition is that for energies E in an interval [a, b] 
and some < s < 1, 

E(\G(x,y;E + ir])\ s ) < C s e^ x ~ yl , (1.6) 

for all rj 7^ 0. Here and henceforth E represents the average over the randomness and 
C < oo and /j, > are constants which may change from line to line, but are to be 
understood as independent of rj. The value of s is of little consequence (if the condition 
(1.6) holds for some s then by Holder inequality it extends to all smaller s > 0), but the 
restriction s < 1 permits to avoid the divergence explained below. 
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E (energy) ► 



Figure 1: Different regions (schematically) in which localization occurs for an operator 
of the form H = H + \V ran dom'- 1) high disorder, 2) extreme energies, 3) weak disorder, 
away from the spectrum of H Q , and 4) band edges. The fractional moment methods 
have been developed for the first three regimes. 



The condition (1.6) was established for a broad class of systems, in any dimension, 
under any of the three conditions: 1) high disorder, 2) extreme energies [4], and 3) weak 
disorder [5] away from the spectrum of the unperturbed operator (A = 0), see Figure 1. 
Localization is known also to occur at the band edges (case 4), for which it can be proven 
[6, 7] by the multiscale approach of Frohlich and Spencer [8]. However, in this more 
delicate situation condition eq. (1.6), which leads to the implications discussed below, 
has not been established yet. Neither has the condition been derived in the continuum 
(for which localization results can be found in refs. [9, 10, 11, 12, 3, 13, 14, 15]). 

We shall recapitulate below why eq. (1.6) can be viewed as a natural technical 
expression of localization, and present a heuristic derivation along the lines of [4]. First, 
however let us list some of its implications. 

l.b Physical implications of the resolvent condition 

There is a growing list of readily identifiable physical properties of an electron gas 
which follow from eq. (1.6), some of which have not been derived without it. 

A new statement which is added here to that list is the exponential decay of the 
two-point function in the ground state |0)) if the Fermi energy, Ep, falls within a range 
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of energies for which eq. (1.6) holds. In terms of Fock-space Fermionic operators 

E(\((0\^(x)4>(y)\0))\)<Ce-^-y\ , (1.7) 

or, in terms of the relevant one-particle spectral projection P<e f , on the energy range 
(-00, E F ], 

E(\{x\P< EF \y)\) < Ce~^ . (1.8) 

Remark: i) If E F falls in a band of extended states the above kernel decays by only a 
power law. 

ii) The derivation of the exponential decay, given below, does not require eq. (1.6) 
to hold for all the energies below Ep . Thus, it applies also to the case in which the 
Fermi level is in a localized regime above a number of bands of extended states. 

iii) The decay presented in eq. (1.8) has interesting implications on electrical con- 
ductance, both in the absence and in the presence of magnetic field (Hall conductance), 
which are discussed below. 

Before we turn to the derivation of the conditions (1.6) and (1.8), let us list some 
other physically meaningful implications of eq. (1.6), some derived by other authors, 
which it may be useful to have listed together. These include: 

• Pure-point spectrum: the spectrum of the operator H in the interval [a, b] almost 
surely consists only of (non-degenerate) eigenvalues with exponentially localized 
eigenfunctions. 

The implication is through either the dynamical localization expressed in eq. (1.9) or the 
Kotani argument [16], as further explained by Simon- Wolff [17]. This argument yields 
the useful principle that those decay properties of the resolvent which hold for almost 
all energies in some interval — in a sense which is not affected by randomizations which 
refresh the site potentials — are typically manifested also by all the eigenfunctions the 
operator may have in that interval. (The property of interest here is the exponential 
decay.) 

• Dynamical localization: assuming eq. (1.6), wave packets of states with energy in 
the range [a, b] do not spread. The key estimate is ([5]) 

E(sup\(x\P [aM e- itH \y)\)<Ce-^ . (1.9) 

One may note that eq. (1.9) is a stronger statement than the exponential localization 
of eigenfunctions (it was not available before eq. (1.6)). Implicit in its derivation is 
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an extension which permits to replace e ltH by an arbitrary bounded function f(H). 
Expressed in terms of the eigenfunctions, with energies E n the assertion is: 

e{ \M*)\\My)\) < ce-^ . (1.10) 

E„e[a,t] 

• Absence of level repulsion: Minami [18] proved that the localization condition 
eq. (1.6) implies that the local distribution of the energy levels in the interval [a, b], 
for the system in [— L,L] d , converges (as L — > oo) to the Poisson law — i.e., the 
energy levels appear to be independently distributed. 

The decay presented in eq. (1.7) has interesting implications on conductivity, both 
in the absence and in the presence of magnetic field. 

We refer here to the conductivities as given by linear-response calculations. We 
shall not address here the interesting questions concerning the validity of such approxi- 
mations, and the role of edge currents. 

• Vanishing of the d.c. electrical conductivity in the absence of magnetic field: We 
shall see below that, for any dimension, 

condition (1.6) =>• a itj {E F ) = (for E F e (a, 6)), (1.11) 

where a is the d.c. electrical conductivity of an electron gas with Fermi energy E F , 
at the zero temperature limit and at zero magnetic field, based a linear response 
calculation (Kubo formula): 

2 

a h3 (E F ) = \im^- ]T x iXj E{\G{^x-E F + iri)\ 2 ) . (1.12) 

Let us remark that this expression for the conductivity follows from the more standard 
Kubo formula ([2, 19]) under the assumption of finite conductivity. (For completeness, 
we present the argument in Appendix A.l.) An earlier proof of the vanishing of a (in 
this form) was provided in ref. [8] for the regime covered by the "multiscale analysis". 

In the presence of magnetic field, one is interested in the Hall conductance. The 
linear response calculation (see Appendix A. 2) is facilitated by the condition 

E( J2 \x\ 2 \(0\P<e\x)\ 2 ) <oo, (1.13) 

Zd 

which is implied by eq. (1.8). The Kubo formula for this situation is 

a id (E) = iTr(P< E [[ Xi ,P< E ], [xj,P< E ]]) . (1.14) 
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• Integral Quantum Hall Effect (2D): Bellissard, van Elst, and Schulz-Baldes [20] 
(BES) proved that if eq. (1.13) holds for a two-dimensional system then at the 
corresponding Fermi energy the Hall conductance ai^Ep) is an integer xe 2 /h and 
it is constant throughout intervals, like [a, b] of eq. (1.9), in which the localization 
length is is uniformly bounded. 

Results leading to this conclusion were also developed by Avron, Seiler and Simon [21] 
(AS 2 ). 

To put this result in a clearer perspective, consider the continuum case of a Landau 
Hamiltonian weakly perturbed by a random potential: 

1 / i B A x\ 2 

H = 2{h- + e ^2^J + XlJ ™ndo m (x) (1.15) 

where U ran dom{.x) could be of the form U ran dom(x) = J2j i]jV(x — Xj), with {xj} randomly 
distributed points and {rjj} independent random coefficients. Had all the results which 
are proven for the lattice Hamiltonians been true in this case, one could deduce that for 
small A the Hall conductance, as a function of the Fermi energy exhibits several plateaux, 
increasing as cr Ha u = (e 2 /h)n. It would indeed be of interest to see an extension to the 
continuum of the localization analysis discussed here. 

The argument of BES [20] is based on a number of sophisticated results of non- 
commutative geometry (in particular, theory developed by A. Connes). To make these 
results more transparent, we include below a direct and simple derivation of the implica- 
tion of eq. (1.8) for the IQHE. The discussion incorporates ideas which were developed 
by AS 2 [21], discussed here in the context of operators with random potentials under 
the localization condition eq. (1.8). 



2. Localization Bounds for the Resolvent 

To convey the key arguments leading to localization bounds let us review the 
derivation of equation eq. (1.6) for the situation with high disorder, or at extreme energy 
(cases 1) and 2) in Fig. 1). Furthermore, let us consider the case where the magnetic 
field B is either zero or constant, K x>y is restricted to the nearest-neighbor pairs (e.g., 
K XjV = the incidence matrix), and U^ er = 0. The equations defining the resolvent (with 
y = and irj absorbed in E) is: (H — E)G(x, 0; E) = S x ^, or: 

(E-XV x )G(x,0;E) = £ K x , x+n G(x + n, 0; E) - 6 X , . (2.1) 

neZ d ,\n\=l 

The solution of this equation does not propagate well (it attenuates, or decays exponen- 
tially) in regions where E — \V X falls out of the spectral range of the hopping operator 
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seen on the right side, i.e., where \E — \V X \ > 2d. If A is large enough (or E is very large), 
most of the lattice will belong to this attenuation set, and one may expect G(x, 0; E) to 
decay exponentially in \x\, indicating exponential localization. 

The big gap in this intuitive argument is that one still needs to address the possible 
tunneling between the sparse sites at which \E — XV X \ < 2d. This gap was closed by the 
"multiscale analysis" of ref. [8], which was developed to handle the technical problems 
caused by possible resonances, manifested through small denominators. An alternative 
approach is to look at suitable moments of G(x, 0; E), with the hope that the averaged 
quantities will be informative enough. Here the small denominator problem shows up 
as follows: for H taken to be a finite matrix (and E real), 

Av(\G(x,0;E)\) = oo (2.2) 

where the average is either over the energy E, integrated over any interval which in- 
tersects the spectrum, or over the values of V x . The latter singularity is explained by 
"rank-one perturbation formulae" (or, alternatively, Cramer's rule), e.g., 

G(x, x; E) = ^ (2.3) 

V ; G(x,x; E)~ l + \V X V ; 

where G(x, x; E) — the value of G(x, x; E) for V x changed to — does not depend on 
V x . In a system in which V x is independent from the values of the potential at other 
sites, there will be rare resonant situations, in which the denominator in eq. (2.3) is 
very small. Despite its rarity, this phenomenon leads to "1/t tail" in the distribution of 
\G(x, x; E)\, as well as of the other matrix elements, and to the divergence of the mean 
values expressed by eq. (2.2). 

Once the nature of the singularity is understood, one may see how to keep it from 
obscuring the picture. The key observation is that this singularity does not cause blowups 
in fractional moments, i.e., averages of \G(x, 0; E)\ s with any < s < 1. Averaging both 
sides of equation (2.1) raised to such a power < s < 1, one may obtain the following 
relation (with the help of a decoupling argument, discussed in [4, 5]) 

cE(\XV x -E\ s ) E(\G(x,0;E)\ s ) < £ E(\G(x + n, 0; E)\ s ) + 5 xfi , (2.4) 

neZ d ,\n\=l 

with c a finite constant depending on the distribution of the random potential. In the 
above relation the effect of the rare resonances is averaged out, and the simple argument 
indicated below eq. (2.1) can be followed in a conclusive way. When 

2d , . 

i= 7wm^W) <l (2 ' 5) 

eq. (2.4) implies that E(\G(x+n, 0; E)\ s ) is a strictly subharmonic function on the lattice 
(its value at x is less that 7(< 1) times its average over neighbors). This readily implies 
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the exponential decay expressed in (1.6), for strong disorder (A large), and at extreme 
energies. This argument leads to the following result. 

Theorem 1 For a random Hamiltonian as in eq. (1-1), with V x independent variables 
having the probability distribution r(v)dv with r(v) bounded, if for some < s < 1 and 
H > 

A>2||r|U(-^- ]T \K , x \ s e^) 1/S , (2.6) 

then 

E(\G(x,y-z)\ s )<Ce-^-y\ , (2.7) 

uniformly in z <E C \ R. 

The complete derivation is given below in Appendix B, where we reproduce in a 
slightly streamlined fashion the argument of ref. [4], and derive the exponential decay 
in a more general setup, in which the probability distribution is not required to have a 
density with respect to dv. This approach can also be applied to other regimes, where the 
resolvent can be studied through other equations, e.g., the relation to the unperturbed 
resolvent operator: G = G Q — G XVG ([5]), yielding bounds which are uniform in the 
two natural cutoffs: finite volume, and imaginary energy shift irj ([22]). 



3. Exponential Decay for the Two— Point Function 

We now turn to the implication of eq. (1.6) for the two-point function, which for 
temperature T > is given by 

(^(x)^(y)) T = E({x\0 T (H - E F )\y)) , (3.1) 

with the Fermi distribution 6 T (u) = (1 + exp(w/T))- 1 . For T = 0: o (u) = I[u < 0]. 

At T > this function always exhibits exponential decay, but at a rate for which 
the general bound — independent of V and E F — vanishes with T: 

\(x\9 T (H - E F )\y)\ < Ce-^ x ^ , (3.2) 

(as follows from Theorem 3 below). We can say more under the localization criterion 
eq. (1.6): 
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Theorem 2 If eq. (1.6) holds at the Fermi energy E F then the two-point function decays 
exponentially fast at the ground state (filled 'Fermi sea'), and also at positive tempera- 
tures — with the correlation length (£) staying uniformly bounded as T — > 0, i.e. 

E(\{x\0r(H - E F )\y)\) < Ce-W* (3.3) 

for all T > 0, and also 

E{\(x\P< EF \y)\)<C^ x - y ^ (3.4) 
(corresponding to the case T = 0). 

Proof: To extract the information from the resolvents, it is convenient to employ the 
contour integral representation. For the projection P<e f we chose the path to consist 
of segments joining — oo — i, E F — i, E F + i, — oo + i. Splitting the integral, we get 

P< Ef = Q 1 (H - Ep) + Q 2 (H - E F ) , (3.5) 

with 

Qi(z) = ±- f d V —, (3.6) 

Z7T 7-1 17? — Z 



1 r° 



2lTl 



-oo 



1 



Q&) = 7^ / du : . (3.7) 



u — 1 — Z U + l 



(due to the randomness, the probability of there being an eigenvalue exactly at the 
energy Ep is zero ([17]), a fact immediate from (1.6)). 

For T > 0, 0p(z) has poles at i7rTx{odd integers}, which grow dense on the 
imaginary axis as T — > 0, with residues equal (— T). Evaluating the Cauchy integral 
along the boundary T of the strip |Imz| < rj with r? = [^f] 27rT 1), where [X] denotes 
the integral part of X, one finds 9 T (H - E F ) = 6 T ,i{H*- Ep) + 9 Tj2 (H - Ep) with 

«tM = T E ~h- z , (3.8) 

nodd;\nirT\<r) 
If 1 

0t,2{z) = — / dwd T (w) 

Z7T1 7r 



27ri 7r w — z 

1 r°° r 1 1 i 
= — / du6 T {u) : : . (3.9) 



2vri 



oc 



U — Yt] — Z U + YX] 



The sum defining 9 T ,i(H — E F ) looks, at small T, like a discrete approximation of the 
integral seen in Qi(H — E F ). 

In each case (T = and T > 0), one may expect the first term to be the more 
delicate one, since it involves resolvents at arbitrarily small distances (of the complex 
energies) from the real axis. Indeed, it is at that point that the assumed condition 
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eq. (1.6) enters. For the more regular terms, Q2(H — E F ) and 9t,2(H — E F ), a starting 
point is provided by the Combes-Thomas estimate (reproduced below in eq. (D.3)), 
however we need to improve on that in order to address the question of the convergence 
of the resulting integrals. (This question can be avoided in case the potential is uniformly 
bounded (H > E > — oo) by closing the contour at any point below E .) 

More explicitly, we estimate the first term as follows: 

EQixlQ^H - E F )\y)\) < ±- £ dr]E(\G(x, y; E F + i v )\) 

< — f 1 dr)r]-^E(\G(x,y; E F + ir])\ s ) < Ce"^"^ , (3.10) 

where we combined the assumed resolvent condition eq. (1.6), with the general operator 
bound \G(x,y;E F + irj)\ < r/" 1 . A similar estimate holds for E(\{x\9 T ,i(H)\y)\), the 
difference being that the integral over rj is replaced by a corresponding Riemann sum. 

The exponential decay for the corresponding second pair of terms, E(\(x\Q2(H — 
E F )\y) |) and E(\{x\9t,2(H — E F )\y)\), is a direct consequence of the following general 
result, whose proof is given here in Appendix D. 

Theorem 3 Let H be as in eq. (1-1) and F be a function analytic and bounded in the 
strip \\vaz\ < rj (by WFW^). Then 

\{x\F{H)\y)\ < W2HFIU e-"!*-"! (3.11) 

for any ji such that the quantity b(fi) = Y< x£ z d \Ko,x\ (e^' — 1) satisfies b(fi) < r//2. 

In fact, this result relies on F having a representation (D.2) analogous to (3.7, 3.9). 

□ 



4. Hall — Kubo conductance as a charge— transport index 

The analysis described above applies in particular to systems with a uniform mag- 
netic field. This case has, of course, attracted a great deal of attention due to the 
remarkable phenomenology associated with the Quantum Hall Effect (QHE). The Inte- 
gral QHE ([23]) (unlike the Fractional case [24]) is understood now to be accountable 
for within the electron-gas picture, in which the particles (or excitations) are subject to 
a one-particle effective Hamiltonian of the type considered here (see, e.g. [25]). 

It has been pointed out that under suitable circumstances the values of anh/e 2 
express topological indices, which would account for both the observed integer values 
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and for the robustness of the phenomenon of IQHE ([26, 27, 28, 29]). Curiously, the ro- 
bustness is re-expressed in the fact that similar conclusions are reached through different 
explanations, in which the topological aspect of Hall conductance appears in different 
disguises. 

In this section we recount one of the approaches to Hall conductance in two di- 
mensions, employing the charge transport index which was introduced and used very 
effectively by Avron, Seiler and Simon (AS 2 ) [21]. The only addition in this paper to the 
above work is the derivation of the exponential decay of the kernel of the projection oper- 
ator P<e f , eq. (3.4), which (in a weaker form, e.g. fast enough power law) is essential for 
the integrality of Hall conductance. Another approach, which was developed by Bellis- 
sard, van Elst, and Schulz-Baldes (BES) [20] will be mentioned in the next section. BES 
proceeded along a slightly different path, employing the Chern-character view of the 
Hall conductance and theorems proven in the context of "non-commutative geometry" . 
The latter work had a more intense focus on random systems, and stated conditions un- 
der which Hall plateaux exist. However, as noted in both works, the seemingly parallel 
tracks actually meet, through a formula discovered by A. Connes. 

The first step may be the formulation of a mathematical expression for the Hall 
conductance within the model considered here. One intriguing option is based on the 
charge pump mechanism proposed by Laughlin [30]. Consider a system in which the 
charges are confined to a plane (e.g., a suitable interface), and the magnetic field is 
changed through an adiabatic process which results in the increase of the flux through 
a finite region D by A$. Changes in the magnetic field are accompanied by an electric 
field (E_), including in the area surrounding D, and current — whose density we denote 
by J_. The rate of charge transport across a contour C encircling D is 



where <7d and an are elements of the bulk (homogenized) conductivity tensor (within 
the plane) 



The last integral on the right side of eq. (4.1) (the induced electromotive force) is tied, by 
Lenz's law, to the flux change —d^/dt. The first term vanishes in situations in which the 
direct conductivity (<7d) is zero. In that case, the integral over time yields an expression 
for the Hall conductance as the ratio of the transported charge to the flux change: 




(4.1) 




(4.2) 



AQ 



(4.3) 



a H = 



A$ ' 
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One may note that it might be easier to analyze increments of flux in multiples of h/e, 
since the addition of such a flux quantum can be accomplished by means of a gauge 
transformation, e.g., 

U a ip(x) = e~ i8a{x) 4>(x) (4.4) 

where a (<£ Z 2 ) is the location of the added flux line and 9 a (x) is the angle of sight from 
a to a; (Arg (x — a), in the terminology of the complex plane). The natural geometry 
for a charge pump based on this mechanism is the Corbino disk, where the transfer 
occurs between two conducting rings, with the region between them filled by material 
whose microscopic structure is modeled by the system discussed in this paper. Detailed 
analysis of this effect was presented in the works of Laughlin [30] and Halperin [31] . 

Motivated by considerations related to the above discussion, Avron, Seiler and 
Simon [21] proposed an interesting representation of the Hall conductance in the model 
discussed here, in the T = limit. They prove that if the two-point function (x\P<E F \y) 
decays fast enough, then in a well defined sense only a finite number of states are moved 
across the Fermi level. The mathematical expression of this is an index, which for a pair 
of projections P and Q of compact difference is defined as: 

Index(P,Q) := dim G H | ^ = ^ } - dim|^G^| £^ = ° j , (4.5) 

(if (P — Q) is a compact operator, the above dimensions are finite). 

Assuming that the above charge transport index coincides also with the charge 
transferred in the course of an adiabatic transition from H to UHU*, (A$ = h/e), AS 2 
take for the Hall conductance the quantity 

a H = j^E(lndex(P< EF ,U a P< EF U* a )) . (4.6) 

The AS 2 study of this quantity rests on the following gem which they added to the 
theory of Hilbert space operators. 

Theorem 4 ( [21]) Let P and Q be a pair of orthogonal projections in a separable Hilbert 
space Ti, whose difference P — Q is a compact operator. If for some integer n > the 
operator (P — Q) 2n+1 is trace class, then (with no further dependence on n) 

tr (P - Q) 2n+1 = Index (P, Q) . (4.7) 



This fact has a simple explanation through the observation that the spectrum of 
P — Q (which consists of a collection of proper eigenvalues in the interval [—1,1]) is 
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symmetric under sign change — except for possible eigenvalues at ±1. A particularly 
elegant proof can be found in [32]. 

Further properties of the index are: 

i) additivity - 

Index(P, Q) + Index(Q, R) = Index(P, R) , (4.8) 
for projections P, Q, R which differ by compact operators, and 

ii) stability - 

Index(P, Q) = Index(P, UQU*) , (4.9) 
under unitaries U with compact difference (U — I). 

(AS 2 prove the above statements by reformulating Index(P, Q) as the Fredholm index 
of PUP in Range P, and invoking known properties of the latter.) 

Two unitary operators U a and Ub, which differ only at the location of the extra 
flux line, are equivalent as far as the Hall conductance eq. (4.6) is concerned, since: 

i) (UaUfr 1 — I) is a compact operator and, by implication, 
ii) 

lndex(U a PU* a ,U b PU;) = 0. (4.10) 

It follows that the charge transport index does not depend (neither in its existence nor 
in its value) on the location of the extra flux line. 

The statement that some power of T = P — U a PU* is trace class may be verified 
by making use of: 

Lemma 1 For an operator with the matrix elements T x ^ y 



|T|| 3 = (tr|T| 3 ) 1 /3< £ J2 \ Tx+b J* 



beZ 



(4.11) 



WZ 2 J 



Proof: One may apply the norm's triangle inequality to the decomposition T = 
J2 be z 2 T^ b \ where = T Xty 5 x ^b,y For each operator in this sum: ||T^||3 = 

|| T (6)* T (fe)||V2 where (r(b)* T (b)} xy = \T x+btX \ 2 5 Xty is a diagonal operator for which the 
norm calculation is elementary. □ 
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The lemma implies 

E(\\Th)< E 2 (E 2 ^(I^ + mI 3 )] • (4-12) 
beZ \xeZ ) 

If the Fermi energy is at a value for which the localization bound eq. (3.4) applies, then 
the condition £?(||r|| 3 ) < oo is satisfied for T XjV = P x>y (l - e~ i(da( - x ' ) ' ea( - y ^), as is easily- 
seen from the bound 



C | 


x-y\ 


3 


1 + 


\x — a 


3 



E{\T x , y \ 6 ) < E(\P x , y \) |1 - e -W*W-»<>(v»\* < A e -^ x - y \ -. ^ . (4.13) 



In this situation, the combination of eq. 4.12, Theorem 4 and some elementary algebra, 
imply that: 

i) the charge-transport index is well defined, 

ii) it is given by 

Index(P<£ F , U a P< EF U* a ) = tr (P< Ep - U a P< EF U* a f (4.14) 
= 2i Yj P x, u P u ,vPv,x [sm(l(u,a,x)) + sin(Z(u, a, u)) + sin(Z(x, a, v))] 
x,u,veZ 

with a an arbitrary point in R 2 \Z 2 , 

iii) the above takes an integer value (which does not depend on a). 

Furthermore, as noted (in slightly different contexts) by Connes [33], BES and 
AS 2 , Index ( P< e f , U a P<E F U*) is a translation invariant function of the randomness (a 
consequence of equations (4.10) and (4.8) ). Since this function is also measurable and 
integrable, Birkhoff 's ergodic theorem implies that the index does not fluctuate, in the 
sense that for almost every realization of the random potential it takes the value given 
by its mean. A significant corollary is that even the mean takes an integer value, i.e. 
the Hall conductance, as represented by eq. (4.6), takes the values 

e 2 

Oh = t~ n 7 with n G Z (4-15) 

and is given by the formula 

e 2 

a H = — 2i E E ( p x,uP u ,vPv,x) (sin a + sin/3 + sin 7) , (4.16) 

x,u,v<=Zj 
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or, using translation invariance, 

e 2 

a H = -r 2i X! E(P 0tU P U!V P Vt0 )(sma + sin/3 + sin 7) , (4.17) 

where 7} are the angles described explicitly in eq. (4.14), and in the second ex- 

pression the summation over x is replaced by a sum over a, varying over the dual lattice. 

The above discussion leads also to the statement, formulated by Bellissard et al. 
[20], that the Hall conductance is constant in regions in which a localization estimate, 
like our eq. (1.8), holds uniformly in E. To see this result, it is convenient to first relate 
the above expression for au with the other expression which was proposed for it, in what 
is known as the Streda formula [34]. 



5. Relation with the Streda-Kubo formula 

In his work on non-commutative geometry, A. Connes [33] presented a remarkable 
formula, whose discrete version reads as follows. For any u,v G Z 2 : 

[sin(Z(u, a, 0)) + sin(Z(i>, a, u)) + sin(Z(0, a, v))] = ir (u 2 vi — u\v 2 ) 

Z2* 
_ 

= Ti u A v . (5.1) 
(For completeness, a streamlined derivation is included here in Appendix F.) 

Using Connes' formula, the expression (4.17) derived for the Hall conductance 
starting from the charge transport index is transformed to (with P = P<e)' 

2m ]T E (P , u Pu,vPv,o) (uAv) (5.2) 

2-ni E ({0\PX 2 PX 1 P\0} - (0\PX 1 PX 2 P\0)) 

2ni E({0\P[[X 2 ,P],[X u P]]P\0)) = a 2>1 (E). (5.3) 

(By translation invariance, the rightmost projection in eq. (5.3) can be omitted.) 
The above expressions are of interest for a number of reasons: 

i) In a suitable setup, the expression provided by eq. (5.3) takes the form of a Chern 
number, and thus it offers another perspective on the topological aspect of Hall 
conductance. (This topic will not be covered here, as it is extensively discussed in 
references [27, 35, 20].) 



e 
~h 



h 
h 
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ii) The above expression coincides with the Kubo formula eq. (1.14) for conductance, 
based on a linear response calculation, e.g., like the one presented here in Ap- 
pendix A). 

iii) The expression provided by eq. (5.2) is very convenient for the derivation of suffi- 
cient conditions for the continuity of the Hall conductance, i.e., for the existence of 
plateaux. The following result is fashioned on a theorem formulated by Bellissard 
et al. [20] (where the assumption is slightly different). 

Theorem 5 (Slightly modified version of a result in [20]) For a random Schrodinger 
operator, H = K x y +XV x , with K incorporating a uniform magnetic field, as in eq. (1.2), 
and V x a random potential whose probability distribution is invariant and ergodic under 
translations, cth(E) (the zero-temperature Hall Conductance at Fermi-energy E) is a 
constant integral multiple ofe 2 /h throughout each interval of energies E, over which for 
some q > 2 the quantity 

6? = E E (\(0\ P<e \x)\ q ) 1/q \x\ (5.4) 

is uniformly bounded. 



Z 2 



This statement depends on the continuity of the integrated density of states, 
-E((0| P<e |0)), a fact which is known for all translation invariant random operators 
in the setting eq. (1.1), [36]. 

Proof: By an elementary telescopic decomposition, and an application of translation- 
invariance, eq. (5.2) implies, for l/q + l/q + l/r— 1 (r > 1), 



cr H (E + AE) - a H (E) | < (5.5) 
2e* 
h 

1/9 



< 3? E «(ii*i«) 1/, *(iiti«) I/, *(i^n I/r mi«-«i 



< 



u,v£Z 2 

6, 2 

I u| 



h 



£(|AP ,o|) 1/r (5.6) 



where P* is either P<£ or P<e+ae, 

AP = P<e+ae — P<e = P(e,e+ae]i and use was 
made of the Cauchy-Schwarz inequality |AP„ i0 | < |APo,o|- The last factor on the right 
side tends to zero by the afore-mentioned general continuity results, while the other 
factors stay bounded under the assumption eq. (5.4) that the localization lengths stays 
uniformly bounded. □ 

An explicit estimate showing the continuity of the integrated density of states 
(though in less than full generality) is the Wegner bound [37]: 

£(|AP ,o|) < A- X |AP| IMU, (5.7) 
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which is valid for random Hamiltonians where the probability distribution for the po- 
tential has a bounded density function r(V) < \\r\\oo. 



Appendix A. The Kubo Formula for the Electric Conductivity 

In this Appendix we present the linear response calculations leading to the expres- 
sions we invoked for the conductance. In particular, we shall reconcile the expression 
(1.12) for o with another familiar form of the "Kubo formula". We shall not address 
here the question of the validity of the linear response theory, which requires a more 
thorough analysis. 

To derive the Kubo linear response formula for conductivity in a system of non- 
interacting Fermi particles consider switching an electric field E_ adiabatically through 
the time-dependent Hamiltonian 

H{t) = H + E-xe vt , (-00 < t < 0, rj | 0) . (A.l) 

The unperturbed density matrix p shall be in equilibrium w.r.t. H, i.e., [p,H] = 0. A 
typical example is the Fermi distribution 

p = 6 T {H-E F ) , (T>0). (A.2) 

The perturbed density matrix p(t) satisfies the initial value problem 

j t p{t) = -i[H(t),p(t)] , t ^J m p{t)e- im = p . (A.3) 

To first order in E (the linear response theory) the solution to eq. (A.3) is 

p(0) - p = -i [° dte^e lHt [E ■ x, p}e~ mt . (A.4) 

For the current density due to the field this yields 

j = -Tr(v(p(p)-p)), (A.5) 
where v = i[H,x\ is the velocity and Tr denotes the trace per unit volume: 

Tr A = limJA^ 1 J2( X \ A \ X ) ■ ( A - 6 ) 

AlZ X £A 

Therefore ji = CijEj, with 

(Tij = - lim Tr( f° dte vt e- iHt [H, x^^xj, p\) . (A.7) 

»?10 J-oo 
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Translation invariance permits one to replace here Tr by an average over the disorder of 
the diagonal term: 

TtA = E((0\AJ0)). (A.8) 

The ergodic argument enabling eq. (A.8) was presented in similar context in 
ref. [20]. Let us consider the probability space whose points are the random environ- 
ments, i.e. potentials u = {V x } ^d, and let T a uo = {V^_ a } be the shift by a e Z d . 
We note that T a act as ergodic shift. An observable A = {A^j^n is stationary if 

U{a)A w U{a)- 1 = A TaW , (A.9) 

for all vectors a which are periods of U^ er . Here U(a) are the magnetic translations 
(1.4). The Hamiltonian (1.1) is stationary in this sense. For stationary operators with 
(x\Au\y) = for \x — y\ large enough and all u G Q the trace (A. 6) almost surely 
takes the value given by the right side of equation (A.8). That expression is valid for 
jjper _ otherwise |o) should be replaced by an average over \x) with x ranging over 
a unit cell. Note that this defines a linear, positive, commutative functional of A. This 
functional then naturally defines a trace class of operators, to which eq. (A.8) extends 
by continuity. 

We now separate the discussion into two cases: zero magnetic field where the 
system is time reversal invariant, and no n- vanishing magnetic field, in which case we are 
interested in the Hall conductance. 



A.l Time reversal invariant systems 

At non-zero temperatures the density matrix is of the form p = f(H) with a 
smooth function / with ||/'||i < oo. In such case, for time reversal invariant systems, 
where (A. 7) is symmetric in {i, j}, eq. (A. 7) can be brought to the form 

a itj = -7rlimTr( / / 5 V (X - ^) l^JM d p v . dP v ) , (A.10) 
vio v J J A — jj, ' 

where 5 v (x) = (ri/ir)(x 2 + rf)~ l . Equation (A. 10) corresponds to a well-familiar form of 
the Kubo formula ([2, 38]). 

We shall now see that this formula yields the expression seen in eq. (1.12): 

2 

&i,j(E) = lim — Yl XiXjE(\G(0,x;E + iri)\ 2 ) . (A.ll) 

VlO IT . 

(tilde added to avoid confusion). The main assuption we shall use is that <7ij(E) is finite 
at all energies. Rougly speaking, this corresponds to a situation where the motion in 
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absence of the electric field is diffusive, at most. On physical grounds one may expect 
that to be the case in the presence of disorder, regardless of localization. Let us remark 
that a related expression for the conductance can be based on the Einsten relation of 
conductance with the diffusion constant ([19]), however we do not use this relation here. 



Theorem 6 Assume that 



sup rf x 2 E(\G(0,x;E + in)\ 2 ) < Const. (A.12) 
v>o „ d 

with a finite constant which applies for all energies. If the limit in eq. (A. 11) exists for 
all E, then 

*ij = -f f'(E)a id (E)dE , (A.13) 

where f(E) = 9 T (E - E F ) and a id = a itj (T) (T > 0). For T -> 0, in lieu of the 
last assumption we require that the limit rj — > in eq. (A. 11) exists for E in some 
neighbourhood of E F , and is continuous there. Under these assumptions, 

Um^ij = °i,j(E F ) . (A. 14) 



Clearly, eq. (A. 14) is the limiting expression for eq. (A.13) as T j 0, where / 
becomes a step function. However a small clarification may be needed, since under the 
weaker assumption made for T = the limit in eq. (A. 10) may exist only for subsequences 
i] n 0. Nevertheless, equation (A. 14) means that the double limit lim-r^o fhn^o is 
unambiguous. 

Proof: We first address the r.h.s. of (A.13): a it j(E) is the limit as i] j of 
^ £ XiXjEdG^x-^ + irj)] 2 ) 

-^E((0\[^ (H-E + irj)' 1 } [ Xj , (H-E- irj)- 1 }^)) 

71 

„2 

— Tr [(H-E + in^v^H -E + ir])-\H -E- i^vAH -E- in)- 1 } 

7T 

= 71 J J 5 V (X - E)8 v (n - E)m i:j (d\, da) . (A.15) 

where m it j(dX,dii) = Ti{dP<\VidP<^Vj) is the conductivity measure ([38]). Under the 
assumption (A.12), the limit 77 J, may be interchanged with the ^-integration seen on 
the right side of eq. (A.13). We thus obtain 

+ J f {E)o-ij(E)dE = lim tt J J g v (X, n)m id (d\ da) (A.16) 
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with 

9r] (X,u) = JdE f G ds[f{E) - f{s\ + (1 - s)u)]5 v (X - £)5„(/i - £) . (A.17) 

The claim eq. (A. 13) is thus equivalent to the assertion that that the right side of 
eq. (A. 16) vanishes. Let us note that rriij(dX,dp) is a finite measure. Using nothing 
more than the smoothness of /, one can show (we skip the analysis here) that: i) g v (X, p) 
is uniformly bounded, and ii) it tends to zero pointwise as r\ \ 0. Thus, (A. 16) vanishes 
by dominated convergence. 

The zero-temperature limit follows by elementary analysis. □ 

A. 2 Systems with decaying two-point function 

Let us return now to the expression eq. (A. 7) for conductance, and discuss it in 
the presence of magnetic field. Our goal is to replace it by a more explicit formula. We 
focus on the zero-temperature limit, where p = P<e f - The following argument, which 
is related to one presented in [20], applies under the (localization) assumption of rapid 
decay of the matrix elements (0\P<e\x)- 

The Hilbert-Schmidt ideal 

T = {A | TtA*A < oo} (A.18) 

is a Hilbert space with inner product (A, B) = Tr A*B. The linear map on T given 
by Ch(A) = [H, A] is self-adjoint. Its resolvent is seen to be 

(C H + Wr 1 ^) = -i f dte vt e- iHt Ae iHt (A.19) 

for i] > 0. Thus, eq. (A. 7) has the appearance of: 

a itj = -i lim Tr((£^ + i^C^Xi) [ Xj , p\) . (A.20) 

In this form, one is tempted to apply the spectral theorem, which implies that 

\im(C H + ir ] y 1 C H (A)=A (A.21) 

for A G (Ker Ch) ± - However Xi is not even in the space T, and hence neither eq. (A.20) 
nor eq. (A.21) applies. In the following argument this difficulty is resolved through the 
replacement of by [[xi, P],P]- 

At zero temperature p = P<e = P is a projection and we have [xj, P] = P[xj, P] 
(1 — P) + (1 — P)[xj, P]P. The substitution of this into eq. (A. 7) amounts, by cyclicity, 
to the substitution of Xj there by the following expression 

(l-P) Xi P + P Xi (l-P) = [[xi,P],P] . (A.22) 
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Unlike x«, the above quantity is in the space T provided 

Tt([x,P]*[x,P}) = e(Y, \ x \ 2 M p <e\x}\ 2 ) < oo. (A.23) 

Zd 

Furthermore, in that case [[xi, P], P] is also in (Ker £h) ± , since for any B e Ker£# we 
have Tr(B*[[xi, P], P]) = Tr([P, B*][x i: P]) = 0. That makes eq. (A.21) applicable, and 
the conclusion is 

a i:j {E) = -iTr([[ Xi ,P],P][ Xj ,P]) = i Tr(P[[x t , P], [ Xj ,P]]) . (A.24) 

Note that a it j(E) is antisymmetric and that, in particular, the longitudinal con- 
ductivity vanishes. In d = 2 it is an integer divided by 2tt by results of [20, 21] and 
reviewed here in Sections 4, 5. If the Hamiltonian (1.1) is time reversal invariant, which 
requires the absence of a magnetic field, the tensor a it j(E) is also symmetric and hence 
vanishes altogether. 



Appendix B. Exponential Decay for the Green Function 

Following is a rigorous derivation of eq. (1.6) for high disorder, along the lines of 
ref. [4] but with somewhat more explicit bounds. As mentioned already, the argument 
can be extended also to other regimes. We allow the probability distribution to be of a 
more general type than considered in Theorem 1. 

Definition Let < r < 1. A r -regular measure q(dv) is one satisfying 

q[v -5,v + 5] < Const. 5 T 

for all v G R, 5 > 0, in which case we let M T (q) be the optimal (smallest) choice for the 
constant Const. 

Such a measure need not have a density dq/dv = r(v). If it does, with r e L P (R) 
for p — (1 — r) _1 , then M T (q) < 2 T ||r|| p . Following is the localization statement. 

Theorem 1' Let < s < r and /i > 0. // 

A > M T (q)^(C^\K ^\^) 1/a , (B.l) 
where C S>T = (2r) _1 (r — s), then 

E(\G(x,y;z)\ s ) < Ce"^^ 1 , (B.2) 

uniformly in z <E C \ R. 
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We begin the proof by stating the following auxiliary fact which plays the role of 
the decoupling lemmas of ref. [4, 5]. Its proof is given in the next Appendix. 



Lemma 2 Let < s < r. Then 

for all r-regular measures dq, ^ dq > and all a, j3 G C . 



Below we will also need a simple upper bound for the r.h.s. We split R into 
\v — [3\~ s < A and its complement. This gives 



J dq{v)\v-(5\~ s <\j dq(v)+J™ d\'q[\v-f3\- s > A'] < -^-M T (q) S ' T (/ dq{% 



l-(s/r) 

(B.4) 



where we minimized over A > 0. 



For the following argument it is important to know that the resolvent is a simple 
rational function of each of the the potential parameters (V x ) at fixed values of the others. 
For matrices that is easily seen from Cramer's formula. More generally, let P x = \x)(x\ 
and let H be the Hamiltonian (1.1) for V x changed to 0. From the resolvent identity 

(H - zY l = (1 + V X (H - z)- l P x ){H - z)- 1 (B.5) 

we get 

nl \ 1 G{x,y;z) 

G{x,y;z) = ~ TT ' ^7 T (R6) 

V x + G (x, x; z) 1 G [x, x; z) 

A simple application thereof is 

sup E(\G(x,x;z)\ s ) < oo . (B.7) 
zeC\R 

In fact already the expectation w.r.t. V x is uniformly bounded by (B.4). 

Proof of Theorem 1': With no loss of generality we set y = and consider eq. (2.1) 
or, more precisely, its replacement for the more general situation (1.1): 

(z - XV X - U*T)G{x, 0; z) = £ K x , y G(y, 0; z) - S x , . (B.8) 

Zd 

To ensure existence we took the resolvent at energies z G C \ R. Raising eq. (B.8) to 
the power < s < 1 yields 

\z-\V x -Ur\ s \G(x,0;z)\ s < £ \K x J'\G(y, 0; z)\ s , (x ? 0) . (B.9) 

Zd 
- 
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Note the particular dependence (B.6) of G(x, 0; z) on V x . Upon taking expectations and 
using eq. (B.3) we obtain 

aE(\G(x,0;z)\ s )< ]T \K x J s E(\G(y,0;z)\ s ), (x ^ 0) . (B.10) 



Zd 

with a = C s>T M T (q)-^ T \ s . 

When a > € g d \Ko,y\ s , the above is a subharmonicity statement for the function 
g(x) = E(\G(x, 0; z)\ s ), which combined with uniform boundedness and exponential 
decay of |-fCr, y | s is known to lead to exponential decay. Following is one of the many 
methods to reach that conclusion (another can be found in ref. [4]). It is based on 
subharmonic comparison. 

For a provisional uniform bound let us note that \\(H — ^r) 1 1 1 < llmz]" 1 yields: 

g(x) < \lmz\- s . (B.ll) 

Thus, g(x) can be viewed as an element in the space of bounded functions £°°(Z d ), and 
eq. (B.10) can be recast as 

a g(x) < (hg)(x) , (x ^ 0) , (B.12) 

where h is the operator with the kernel |i£a; )3/ | s . Note that if y? G i°°{Z d ) with y?(0) < 
satisfies 

cup(x) < (ktp)(x) , (x^O) (B.13) 

with a > C = sup^g^d J2 7 d \K x>y \ 8 , then (fi(x) < 0. In fact, if M = sup yd <p(x) > 
then a<f(x) < CM and we get a contradiction by taking the supremum over x. We apply 
this conclusion to (p(x) = g(x) — g(0)e~^. The length scale /i -1 is set by the condition 

J2 \K ^\ < a , (B.14) 

Zd 

which is eq. (B.l). Using 

(he-M)(x) < ( £ \K J s e^)e~^ (B.15) 

Zd 

we see that aip(x) < (hip)(x) for x ^ and hence (p(x) < 0, i.e., 

g(x) < ^(O)e"^ 1 . (B.16) 
The claim follows now by combining this with eq. (B.7). □ 
For certain applications the following variant of eq. (B.2) is useful: 
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Corollary 1 Under the assumptions of Theorem 1 we have 

G(x,y;z) 



G(x, x] z) 



') < e -^ x -y\ , 



(B.17) 



for all z e C \ R. 



Proof: This is actually a corollary of the proof of Theorem 1. Due to G(x,y;z) = 
G(y,x;z) we may, upon interchanging x and y, prove (B.17) with G(y,y;z) in the 
denominator. We then set y = as before. A short computation based on eqs. (B.5, 
B.6) shows the following dependence 



G(x, 0; z) a 
G(0,0;z) +i5 = V x -/3 



(B.18) 



on V x . If Imz > 0, as we may assume without loss, the regularization by 5 > ensures 
that 

\G(0,0;z) + i5\ 2 = \G(0, 0; z)\ 2 + 25ImG(0, 0; z) + 5 2 > 5 2 , (B.19) 

because of ImG(0,0;z) = Imz {0\(H - z)-\H - ^) _1 |0) > 0. We then divide eq. (B.9) 
by |G(0, 0; z) + i5\ s and obtain eq. (B.10) once more but now for 
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(x) = E( 



G(x,0;z) 



G(0,0;z) + i5 



(B.20) 



This is bounded in x by (5|Im2;|) s . The upshot is again eq. (B.16) with g(0) < 1. 
Hence 



G(x,0;z) 



-fi\x-y\ 



G(0,0;z) + i5 

and the conclusion is by monotone convergence in the limit 5 10. 



(B.21) 

□ 



Appendix C. Proof of the Decoupling Lemma 

We shall need the inequality 



\v-p\~ s + \u-f3\~ s < 



\v-P\< 



\u\- s + \u - (3\~ s ) + 



\u\ 



u-P\ 



-X\v\-' + \v-0\-') (C.l) 



for all u,v, (3 e C (except for vanishing denominators). Multiplication by \v — /3\ s \u — f3\' 
shows it to be equivalent to 



(Si " i) l« " W + H s + + (& - 1) I- - ^ > • 



(C.2) 
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Since this expression is symmetric in u and v it suffices to prove this for \u — f3\ > \v — /3\. 
The triangle inequality yields \u — f3\ s < \v — {3\ s + \u\ s + \v\ s , which we apply to the two 
middle terms of (C.2) so as to bound it from below by 



Is 

— \u 



- W + (ft - 2 ) \ v - (ft + ft " 2 ) \ v - ' 

\u\° x \v\ s J y \u\ s \v\ s J 

f-1 



since t + t > 2 for t > 0. This proves (C.I). Replace there i> by f — a and similarly 
for u, [3, and integrate w.r.t. dq(u)dq(v). The result is 

| dq(u) j dq(v) * < J dq(v) ^ ~ ^ | - + |u - /3|- s ) 

where, actually, each side comes duplicated with dummy variables u, v interchanged. 
The last integral is estimated by (B.4). □ 



Appendix D. Analyticity and Exponential Decay (proof of The- 
orem 3) 

In section 3 we claimed and used the following statement. 

Theorem 3 Let H be as in eq. (1.1) and F be a function analytic and bounded in the 
strip \lmz\ < rj (by \\F\\oo). Then 

\{x\F(H)\y)\ < W2HFIU e^*""! (D.l) 

for any ji such that the quantity b(fj) = Yl, 7* |-^o,x| (e^ — 1) satisfies b(fx) < rj/2. 



By continuity it suffices to prove (D.l) in any smaller strip. We may thus assume 
F to be continuous up to the boundary. We note that under the above assumptions F 
has the representation 



1 



1 1 



F(E) = — du : — f(u-E) 

Z7T1 J l U — IT] U + lT] j 

= D*f(E) (D.2) 

with D{u) = i]/[(u 2 + r] 2 )7i] and / a uniformly bounded function (||/||oo < 00). In fact, 
(D.2) is solved by / = F + + F_ - D * F, where F±(u) = F(u ±ir]T i0). This follows 
from D * (F + + F_) = (D+ + D_) * F and D + + D_ = 5 + D * D. 

The proof of Theorem 3 is related to the Combes-Thomas bound [39] : 

\G(x,y;E + i V )\ < {2/^'^ (D.3) 

with ji as above. In order to integrate over u in eq. (D.2) down to —00 we first develop 
the following related estimate. 
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Lemma 3 With \i be small enough so that b(/i) < i]/2, 

\G(x,y;E + ir))-G{x,y;E-ir})\ < 12 ^e"^-^ 

•(^i (^^ k) 1/2 (yl ^^ iy) 1/2 • (D.4) 

Proof: We set E = for notational simplicity. Let 27rLD = (H — ir/)^ 1 — (H + 177) —1 , 
and for any bounded function f(x) let Df = e f De~ f + e~ f De f . Since 

( e /W-/fe) + e - {f{x) - f{y)) ){G{x, y; irj) - G{x, y; -irj)) = 2m(x\D f \y) , (D.5) 

the desired bound would follow from the statement that for any / satisfying \f(x) — 
f(y)\ < fJ>\ x — y\ ( e -g-; a function which in a suitable finite region is f{u) = ji\u — y\) 

||[#2 + v 2 /2] 1 / 2 D f [H 2 + V 2 m 1/2 \\ < 677/tt . (D.6) 

To prove eq. (D.6), we first group the terms as follows, 

^1111 



Hf — ii] Hf + it] Hj + HJ — ii] 

1111 

-B + B- , (D.7) 



-17] Hf + 17] Hf + 17] Hf - 17] 

with Hf = He~f , B = Hf — H, and 5 ± = (5* — ± 2i^). By the assumption on /i, 
we have 

||-B|| < b{ix) < 77/2 , —377 < i.B ± < 3?7 . (D.8) 

We now claim that 



1 

2 

Indeed, using the positivity of the last term in 



HJ - irj) (H f + irj) > - [H 2 + 7] 2 /2] . (D.9) 



(H* f -irj)(H f + irj) = ^[(H -irj)(H + irj) -2B*B] + ^(H -ir) + 2B*)(H+ir) + 2B) (D.10) 

and eq. (D.8), we see that the l.h.s. is bounded below by (1/2) [i/ 2 + 7] 2 — (t] 2 /2)]. 

The estimates (D.8), (D.9), and equation (D.7), readily imply (D.6). This bound, 
combined with an application of the Cauchy-Schwarz inequality to the r.h.s. of eq. (D.5) 
proves the claim made in eq. (D.4). (The exponential weight is incorporated in the l.h.s. 
in eq. (D.5).) 

Proof of Theorem 3: Let us note that unlike the corresponding integral of operator 
norms, the following integral is bounded: 

du ^ (H-uy + vW X) - <D - U) 
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(using the spectral measure representation). The claim, eq. (D.l), is obtained by com- 
bining the integral representation of F, eq. (D.2), with the exponential bound eq. (D.4), 
and employing the Cauchy-Schwarz inequality to reduce the resulting integral to the 
one estimated in eq. (D.ll). □ 



Appendix E. The ir) Regularization 

The addition of a small imaginary term ir] to the energy is a standard regular- 
ization, and a convenient alternative to the finite-volume cutoff. Such a cutoff appears 
also in the Kubo formula for the electrical conductivity in the absence of magnetic field, 
eq. (1.12). Dealing with such expressions one should bear in mind the operator bound: 
r]\G(0,x;E + irj)\ < 1. For the conductivity given by eq. (1.12) it implies that quite 
generally 

(JiAE) < liminf — V \xixAE(\G(0,x; E + ir])\ s ) , (E.l) 
^° 7T d 

for any s < 2. Thus, the fractional moment localization estimate, eq. (1.6), directly 
implies the vanishing of the Kubo conductivity — in the absence of magnetic field. 

Working with this regularization, it is useful to have also the following lemma. Its 
second bound can be used for yet another proof of the dynamical localization (1.10), 
which was originally derived using the finite- volume cutoff (ref. [5]) and which was 
provided another derivation in ref. [40]. 



Lemma 4 // (1.6) holds, and the probability distribution, q(dv) = r{v)dv, satisfies the 

regularity condition r G L X {R) H L P (R) for some p > 1, then for any E G [a, b] and 

ri^^EUGfay-tE + irfil 2 ) < Ce^ 1 ^ 1 , (E.2) 

/ dEr]E(\G(x,y;E + ir])\ 2 ) < Ce~^ x - y K (E.3) 



To calibrate these statements we note that without any assumptions on the self- 
adjoint operator H: 

/oo 
dEr)\G(x,y] E + irj) | 2 < vr . (E.4) 
-oo 

We shall only sketch here the proof of Lemma 4, which is by arguments seen in 
ref. [5] (Lemma 3.1) and in ref. [22] (Lemma 3). Some of the key points in the analysis 
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are: 

i) Quite generally, for any < s < 2: 

\Imz\\G(x,y;z)\ 2 < \ImG(x, x; z)\ ■ ^ V \% ■ (E.5) 

\(j[X, X, Z) I 

(For s = the proof is by a judicious use the Cauchy-Schwarz inequality, for s = 2 it 
follows from |Imz| < \lmG(x,x; z)~ 1 \, and for other < s < 2 it holds by interpolation.) 

ii) Using eq. (B.6) on the first factor on the right side of (E.5) yields 



|Imz||G(a;,y;*)| 2 < 



\lmG(x,x; z) x | \G(x,y;z) 



\V x + G(x,x;z)- l \ 2 \G(x,x;z)\* ' 

where, again by (B.6), the second quotient is independent of V x (!). 

To derive Lemma 4, one may now average first over V x — in effect making use of 
the high degree of independence of the values of the potential at different sites — and 
then use (B.17). The proof is most direct for the case of bounded density r(v) (p = oo), 
and the extension to more singular distributions is by arguments similar to those found 
in ref. [5]. 

Remark: The condition expressed by the second statement in Lemma 4 implies directly 
the exponential decay law 

E{ sup |(x|P M /(//)|y)|)<Ce-^l, (E.6) 

/:||/l|oc<l 

which is equivalent to eq. (1.10). For that, one may use the resolvents for an approximate 
8 function, writing (for / continuous): P^^f^H) = s—lim^^f^H), where 

f (H) = - [ b dEri f(E) . (E.7) 

JvK ' nJa 'H-E-ir]H-E + ir} Jy ' y ' 

The matrix elements of (E.7), can be easily brought to a form in which (E.3) implies 
(1.10). The key tool is the Cauchy-Schwarz inequality, applied in the different setups: 
the state Hilbert-space, E — the average over the disorder, and in J a b dE. 



Appendix F. Connes' area formula. 

In Section 5 an identity is stated relating two expression for the Hall conductance: 
one based on the charge-transport Index, and the other corresponding to the Streda 
formula which takes the form of a Chern number, eq. (5.3). Following is the derivation 
of Connes' area formula [33] which has been used to prove that relation. The formulation 
and derivation presented here incorporate a streamlined argument of Colin de Verdiere 
([41, 32]), shown to us by R. Seiler. 
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Theorem 7 For a fixed triplet u^ 2 \ G Z 2 , let cti(a) G (— 7r, 7r) fre £/ie ang/e of 
view from a G Z 2 * ofu^ l+2 ^ relative to (with ctj(a) = i/a iies between them). Let 

g{a) be an antisymmetric bounded function satisfying: 



near a = 0. Then, 



g(a) = a + 0(a 3 ) (F.l) 



]T ^^(a J (a)) = 27rArea(A( M ( 1 ), M ( 2 ), M ( 3 ))) (F.2) 

ff* i=l 



where Area (A(. . .)) £/ie triangle's oriented area. 

Of special interest to us is the case with g(a) = sin a, which is used here in eq. (5.3). 

Proof: We may assume the triangle to be positively oriented. The statement (F.2) is 
true for g(a) = a. Indeed, for each a G Z 2 * 



1 

Y,on{a) = 2<nl 1/2 
*=i 



inside 

> for a { on the boundary of 
outside 



> the triangle (F-3) 



Thus, for g(a) = a the l.h.s. of (F.2) is 2nx the number of dual lattice sites within 
the triangle (counting a boundary site with weight 1/2). This number is the same for 
triangles obtained by the lattice translation and reflection symmetry operations. Since 
this set of triangles tiles the plane, the number of enclosed dual sites must equal the 
triangle's 

The above observation reduces eq. (F.2) to the statement that for f(a) = g(a) — a 

E E /(<*(«)) = o • (F.4) 

aeZ 2 1=1 

A significant difference between / and g is that the individual terms /(aj(a)) are 
summable in a G Z 2 , since by eq. (F.l) /(«j(a)) = 0(|a| -3 ) for \a\ — > oo. However, 
each of the three individual sums changes sign under the reflection with respect to the 
midpoint of the corresponding edge, (a^ +1 ^ + a^ +2 ^)/2 G (Z/2) 2 (which is a symmetry 
of the lattice Z 2 ). Thus even the individual sums (at given i) vanish. □ 
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